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Abstract 

Exact solutions of the string equations of motion and constraints are systematically 
constructed in de Sitter spacetime using the dressing method of soliton theory. The string 
dynamics in de Sitter spacetime is integrable due to the associated hnear system. We start 
from an exact string solution g^Q) and the associated solution of the linear system ^('^^(A), 
and we construct a new solution 'J'(A) differing from ^'('^^(A) by a rational matrix in A with 
at least four poles Aq, I/Aq, Ag, l/Ag. The periodicity condition for closed strings restrict 
Aq to discrete values expressed in terms of Pythagorean numbers. Here we explicitly 
construct solutions depending on (2 + l)-spacetime coordinates, two arbitrary complex 
numbers (the 'polarization vector') and two integers (n,m) which determine the string 
windings in the space. The solutions are depicted in the hyperboloid coordinates q and in 
comoving coordinates with the cosmic time T . Despite of the fact that we have a single 
world sheet, our solutions describe multiple (here five) different and independent strings; 
the world sheet time r turns to be a multivalued function of T. (This has no analogue in 
flat spacetime). One string is stable (its proper size tends to a constant for T oo, and 
its comoving size contracts); the other strings are unstable (their proper sizes blow up for 
T ^ oo, while their comoving sizes tend to constants). These solutions (even the stable 
strings) do not oscillate in time. The interpretation of these solutions and their dynamics 
in terms of the sinh-Gordon model is particularly enlighting. 
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1 Introduction 



Since the propagation of strings in curved spacetimes started to be systematically investi- 
gated, a variety of new physical phenomena appeared These results are relevant both 
for fundamental (quantum) strings and for cosmic strings which behave essentially in a classical 
way. 

String propagation has been investigated in non-linear gravitational plane waves and 
shock- waves 0], black holes 0, conical spacetimes @, and cosmological spacetimes [0,0- 

Among the cosmological backgrounds, de Sitter spacetime occupies a special place. This 
is, in one hand relevant for inflation and on the other hand string propagation turns to be 
specially interesting there String unstability, in the sense that the string proper length 

grows indefinitely is particularly present in de Sitter. The string dynamics in de Sitter universe 



is described by a generalized sinh-Gordon model with a potential unbounded from below |14 



The sinh-Gordon function a{a,T) having a clear physical meaning : H~^e°'^"''^^/'^ determines 
the string proper length. Moreover the classical string equations of motion (plus the string 
constraints) turn to be integrable in de Sitter universe |T^,[|1^]. More precisely, they are 
equivalent to a non-linear sigma model on the grassmannian S0{D,1)/0{D) with periodic 
boundary conditions (for closed strings). This sigma model has an associated linear system |^ 
and using it, one can show the presence of an infinite number of conserved quantities In 
addition, the string constraints imply a zero energy-momentum tensor and these constraints 
are compatible with the integrability. 

The so-called dressing method in soliton theory allows to construct solutions of non- 
linear classically integrable models using the associated linear system. In the present paper 
we systematically construct string solutions in three dimensional de Sitter spacetime. We 
start from a given exactly known solution of the string equations of motion and constraints 
in de Sitter and then we "dress" it. The string solutions reported here indeed apply to 
cosmic strings in de Sitter spacetime as well. The dynamics of cosmic strings in expanding 
universes has been studied in the literature for the Friedman-Robertson- Walker (FRW) cases 



(see for example ||^, [T3I4I1H)- must be noticed that the string behaviour we found here 
in de Sitter universe is essentially different from the standard FRW where the expansion 
factor R{T) is a positive power of the cosmic time T. In such FRW universes, strings always 
oscillate in time, the comoving spatial string coordinates contract and the proper string size 



stays constant asymptotically for T ^ cxo [0],||T6[ . In the cosmic string literature this is known 
as 'string stretching'. We called such behaviour 'stable' 0,[0,|jl6[. On the contrary, in de 
Sitter spacetime, as we show below, two types of asymptotic behaviors are present : (i) the 
proper string size and energy grow with the expansion factor ('unstable' behaviour) or (ii) they 
tend to constant values ('stable' strings). 

The unstable string solutions in de Sitter universe may provide a mechanism to self-sustain 
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inflation as proposed in refs.[|TT|-[0] without advocating an inflaton field. The multi-string 
exact solutions in de Sitter spacetime presented here should provide essential clues about the 
feasability of inflationary string scenarios. 

We apply here the dressing method as follows. We start from the exact ring-shaped string 
solution g(o) [0 and we find the explicit solution \Ef*^°)(A) of the associated linear system, where 
A stands for the spectral parameter. Then, we propose a new solution \E'(A) that differs from 
'^^^^ (A) by a matrix rational in A. Notice that ^(A = 0) provides in general a new string solution. 
We then show that this rational matrix must have at least four poles, Aq, I/Aq, Aq, I/Aq, as a 
consequence of the symmetries of the problem. The residues of these poles are shown to be 
one-dimensional projectors. We then prove that these projectors are formed by vectors which 
can all be expressed in terms of an arbitrary complex constant vector |xo) and the complex 
parameter Aq. This result holds for arbitrary starting solutions g(o). 

Since we consider closed strings, we impose a 27r-periodicity on the string variable a . This 
restricts Aq to take discrete values that we succeed to express in terms of Pythagorean numbers. 
In summary, our solutions depend on two arbitrary complex numbers contained in |xo) and two 
integers n and m . The counting of degrees of freedom is analogous to 2 + 1 Minkowski spacetime 
except that left and right modes are here mixed up in a non-linear and precise way. 

The vector |xo) somehow indicates the polarization of the string. The integers (n, m) deter- 
mine the string winding. They fix the way in which the string winds around the origin in the 
spatial dimensions (here ) . Our starting solution g(o) (cr, r) is a stable string winded + rn? 
times around the origin in de Sitter space. 

The matrix multiplications involved in the computation of the final solution were done 
with the help of the computer program of symbolic calculation "Mathematica" . The resulting 
solution q{(J,T) = {q^ , , q'^ , q^) is a complicated combination of trigonometric functions of a 
and hyperbolic functions of r. That is, these string solitonic solutions do not oscillate in time. 
This is a typical feature of string unstability [if- 0- fill- The new feature here is that strings 
(even stable solutions) do not oscillate neither for r 0, nor for r ±oo. Figs. 3-4 depict 
spatial projections {q^,q'^,q^) of the solutions for two given polarizations |xo) and different 
windings {m,n). 

We plot in figs. 5-11 the solutions for significative values of |xo) and {m,n) in terms of the 
comoving coordinates (T, 

T = hog{g' + q') , ^' = ^^ > X^ = ^^ (1) 
H H q^ + q^ H q^ + q^ 

The first feature to point out is that our solitonic solutions describe multiple (here five or 
three) strings, as it can be seen from the fact that for a given time T we find several different 
values for r. That is, r is a multivalued function of T for any fixed a (fig. 5-6). Each branch 
of r as a function of T corresponds to a different string. This is a entirely new feature for 
strings in curved spacetime, with no analogue in fiat spacetime where the time coordinate can 
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always be chosen proportional to r. In flat spacetime, multiple string solutions are described 
by multiple world-sheets. Here, we have a single world-sheet describing several independent 
and simultaneous strings as a consequence of the coupling with the spacetime geometry. Notice 
that we consider free strings. (Interactions among the strings as splitting or merging are not 
considered). Five is the generic number of strings in our dressed solutions. The value five can 
be related to the fact that we are dressing a one-string solution (g(o)) with four poles. Each 
pole adds here an unstable string. 

In order to describe the real physical evolution, we eliminated numerically r = r(cr, T) from 
the solution and expressed the spatial comoving coordinates and in terms of T and a. 

We plot r(cr, T) as a function of a for different fixed values of T in fig. 7-8. It is a sinusoidal- 
type function. Besides the customary closed string period 2tt, another period appears which 
varies on r. For small t , r = T{a, T) has a convoluted shape while for larger r (here r < 5), it 
becomes a regular sinusoid. These behaviours reflect very clearly in the evolution of the spatial 
coordinates and shape of the string. 

The evolution of the five (and three) strings simultaneously described by our solution as a 
function of T, for positive T is shown in figs. 9-11. One string is stable (the 5th one). The 
other four are unstable. For the stable string, (X^,X^) contracts in time precisely as e~^^, 
thus keeping the proper amplitude (e^^X^, e^-^X^) and proper size constant. For this stable 
string (X^,X^) < jj. {l/H = the horizon radius). For the other (unstable) strings, (X^,X^) 
become very fast constant in time, the proper size expanding as the universe itself like e^^ 
. For these strings (X^,X^) > jj. These exact solutions display remarkably the asymptotic 
string behaviour found in refs.0,0]. 

In terms of the sinh-Gordon description, this means that for the strings outside the horizon 
the sinh-Gordon function a{a, r) is the same as the cosmic time T up to a function of a. More 
precisely, 

a{a, t) ^ 2HT{a, r) + log {2H^ [{A\a)'f + {A\ayf] } + 0{e-^"^). (2) 

Here A^{a) and A^{(t) are the X^ and X^ coordinates outside the horizon. For T +oo 
these strings are at the absolute minimun a = +oo of the sinh-Gordon potential with infinite 
size. The string inside the horizon (stable string) corresponds to the maximun of the potential, 
a = 0. a = is the only value in which the string can stay without being pushed down by 
the potential to a = ±oo and this also explains why only one stable string appears (is not 
possible to put more than one string at the maximun of the potential without falling down). 
These features are generically exhibited by our one-soliton multistring solutions, independently 
of the particular initial state of the string (fixed by > and {n,m)). For particular values of 

>, the solution describes three strings, with symmetric shapes from T = 0, for instance like 
a rosette or a circle with festoons (fig. 9-11). 

The string solutions presented here trivially embedd on D-dimensional de Sitter spacetime 
{D > 3). It must be noticed that they exhibit the essential physics of strings in D-dimensional 



4 



de Sitter universe. Moreover, the construction method used here works in any number of 
dimensions. 

This paper is organised as follows: in section 2 we describe the string equations in de Sitter 
universe and its associated linear system. Section 3 deals with the dressing method in soliton 
theory, its application to this string problem and the systematic construction of solutions. In 
section 4 we explicitly describe the starting background solution g(o) (a, r) and the solution 
\E'*^*'^(A) of the associated linear system. In section 5, we analyze our multistring solutions and 
describe their physical properties. 



2 The string equations and their associated linear sys- 
tem 



The string equations of motion in D-dimensional de Sitter space-time can be written in the 
following form: 

Q^v) +l)i(li;Jlv) = (3) 
where q) is a {D + 1) -dimensional real vector of unit pseudolength 

{qJq) = l, (J = diag(-l,l,...,l)) (4) 

and ^,1] are light cone coordinates in the world sheet: 

r = V + t(^ = V-^ (5) 
In addition, we have the string constraints (conformal conditions) 

{q^Jq^) = 0, (g^Jg^) = 0. (6) 

The solution q) should be a periodic function of a = ?7 — ^, with period 27r for closed strings. 
We are going to find solutions of this equation by using the Riemann Transform Method 
|13|. The most important observation is that equation can be rewritten in the form of a 



chiral field model on the Grassmanian Gd = SO{D, 1) /0(D). Indeed, any element g G Go can 
be parametrized with a real vector q) of the unit pseudolength 

g = l-2g)(gj, {qJq) = l. (7) 

In terms of g, the string equations (0)-(§) have the following form 

2 gc,7 + g? g + g % = ' (8) 
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and the conformal constraints are 



trg| = 0, trg2 = 0, (9) 

which are equivalent to eqs.(^). The fact that g G G/? imphes that g is a real matrix with the 
following properties: 

g = Jgh, g2 = J, trg = 2 geSL{D + l,R). (10) 

These conditions are equivalent to the existence of the representation (0). 

Equation (j^) is the compatibility condition for the following overdetermined linear system: 

(11) 



where 



Or in terms of vector q) 



U = gig, V = gr,g . (12) 



U = 2q^){qJ-2q){q^J, 
V = 2q^){qJ-2q){q,J. 
In order to fix the freedom in the definition of \E' we shall identify 

^(A = 0) = g. (13) 

This condition is compatible with the above equations since the matrix function \E' at the 
point A = satisfies the same equations as g. Thus the problem of constructing exact solutions 
of the string equations is reduced to finding compatible solutions of the linear equations (p!T| ) 
such that g = \I/(A = 0) satisfies the constraints eqs.(|^) and (p!0D . 

3 The Dressing Method in Sohton Theory 

3.1 The reduction group of the associated Unear system 



We will consider now the symmetry group (or the so called "reduction group" 0, |jT3[) 
enjoyed by the linear system of equations 



= *, = ^*. (14) 



when eqs.(|TOD hold. 
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It follows from the condition {qJq) = 1 that the matrix g = / — 2q){gJ anticommutes with 
U and V: 

gU + Ug = 0, gV + Vg = 0. 
This implies that the matrix function g\l'(l/A) satisfies the same equation as ^^(A) : 

[^?^(1/A)]c = Y^[^7^(l/A)], [9^il/\)]r, = j^[9^{l/\)l (15) 

Then, it can differ from \I^(A) only on a matrix multiplier which does not depend on ^,,1]: 

g^(l/A) = ^(A)5i(A). (16) 

The vector q), the corresponding matrix g and the currents U, V are real. Therefore \E'*(A*) 
is a solution of equations (|TT]) as well, and we have 

^*(A*) = ^(A)52(A) (17) 

In addition, by using eq.(|l3) twice, we find 

52(A) 6;iX*) = I (18) 

The fact that g G S0{3, 1) yields JU^ J = -U, JV^J = -V and implies that (J^*(A) 7)-^ 
obeys the same equation (|TTD as \I'(A) 

{J^\\)J)-' = ^i\)Ssi\). (19) 

The transformations (|T6|),(|T^) and (|1^) generate a finite group which is called the reduction 
group of the problem and which guarantees that the properties (0) hold for g = \I'(A = 0). 

3.2 Rational dressing 

Suppose we know a particular solution g(o)(^, of ^^e string equations (H). We shall denote 
by t^(o)(^) 0) ^(o)(^; its corresponding currents (p!2D, and by \1/^°''(A, r], ^) , the corresponding 
compatible solution of the overdetermined system ([TTl) . We assume that as well as ?7(o) 
and V(o) are explicitly known. 

To construct a new solution g we assume that the corresponding function differs from 
on a rational matrix multiplier ^(X,r],^) 

^(A) = <I>(A)^(°)(A). (20) 

We assume that $(A) is rational in A, but, of course it might have a complex dependence in ^, 77. 
The dressing method consist on finding a matrix $(A) such that ^I/(A) given by eqs.(|20|) fulfils 
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the linear system (|T4D and the symmetry conditions (p!6|)-([T9|). Then, once $(A) is known, the 
string solution g{ri,^) follows from eq.([T3|). 



It follows from and ( ^01) that $ should obey the following symmetries: 

$(0)^(°)(0)<I>(1/A) = <l>(A)^(°)(0), (21) 
$*(A*) = $(A), (22) 
j$t(A)J = $-i(A) (23) 

We assume that the constant matrices 5i(A), 52(A) and S^^X) coincide for the dressed and 
the undressed solutions. 

Suppose that the rational function $(A) has a pole at the point Aq. It follows from eqs. (pT]) 
and (|2^ that it must have poles at the points I/Aq, Aq, I/Aq as well and, in addition, $(00) = /. 
Thus, the simplest (generic) possible case is 

A A* B B* 

A — Aq a — Ag A — Aq ^ A — Aq ^ 

where A, B are matrix functions of (^, rf) to be determined below. This simplest case will be 
called the one-soliton solution from now on. We choose this name since in the context of non- 
linear integrable equations in an infinite space interval (the sine-Gordon equation, for instance), 
this minimal pole-structure (the minimal number of poles compatible with the symmetry group) 
in A generates the one soliton solution (see for example ||T7| ). 

Here we have taken into account eqs.(^) and (|21|). It follows from eq. ( pSj) that 

_ ^ ^ JAh ^ JAU ^ JB^J ^ JB^J ^25) 

A — Aq a — Aq a — Aq ^ A — Aq ^ 

here f denotes Hermitian conjugation of a matrix. The condition Vt{\) = $(A)$^^(A) = I can 
be imposed in the following way: the right hand side (f2(A)) is a rational function of A which 
takes the value I at the point A = 00, then Q{X) will be identically / if it does not have any 
singularity on the Riemann sphere of A. Double poles would vanish if and only if 

AJA^ = 0, BJB^ = 0. (26) 

Thus the matrices A, B are degenerated and we can write them as a sum of bivectors 

A = J2a,) (x, J, B = J2b^) {y^J. (27) 

i i 

The constraints (EBI) imply 



{xi\J\xj) = for all pairs i,j. (28) 

which means that the vectors Xj) are null and mutually pseudo-orthogonal. Therefore, since 
pseudo-orthogonal null vectors are proportional, we have 
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and without loss of generality, we take: 



A = a){xJ, B = b){yJ. (29) 

Now the constraints (^) read 

{xJx) = 0, {yJy) = (30) 

In addition, by requiring the residues of f2(A) to vanish at the points Aq, I/Aq, Aq, I/Aq , 
we get 

^ T T A ^ , , B J A* J B* J , . ^ 



BJ + JB^ + BJ{-- — T + — ^ — T + 



At _ AJ B*J ^*J_^Bt^Q 



-^0 — Aq ^ Ao — Aq ^ Aq ^ — Aq ^ Aq — Ag 

. A^ B^ , , B*J AJ BJ 

A*J+JA^ + A*Ji r-^ +T^^ + T^^) + ( w A-^ + A^^ + V^^^^ = ° 

^0 ~ ^0 ^0 ~ ^0 — Aq Aq — Aq Aq — Aq Aq — Aq 

/it Rt /it 

B*J + Ji?t + i?V(--^— ^ + + - 



\* \*— 1 \*— 1 \— 1 \* \— 1- 

Aq ~ Aq Aq — Aq Aq — Aq 

A* J BJ AJ ,_ 



Aq ~ Aq Aq — Aq Aq — Aq 

Later on we shall demonstrate that the periodicity condition on cr can be satisfied only in 
the case where all poles (Aq, Aq, Aq ^, Aq~^) of $(A) are purely imaginary [see eqs.(|59D,(|68D and 



(ffOl)]. From now on, we shall denote Aq = in, k E R. Substituting the bivectoral representation 
( P^D in the above equations and by separating bivectors we get the system of vector equations 
with respect to a), b) 

^^.i{x*Jy*) .i{x*Jx) ,,i{x*Jy) ^. 
b)- ^ + — - + b)- ^ = -X*) 

Kj ~\~ Kj 2iK K K 

_a)^^-V)'-^ + a*)'-^^ = y) 

K + K ^ 2k ^ K — K ^ 

By bivectoral separations we mean the following trick: suppose ve have an equation of the form 
p) (X + X) (p = where p),X) are vectors and p) ^ 0. Then the only solution of this equation 
is X) = and the matrix equation is reduced to a vector one. 
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One can solve this system of linear equations and express the vectors a), b) in terms of x),y) 

a) = ^""^""l ~ [{^' - l)x*){y*Jy) + 2(1 + ^')y*){x*Jy) + 2(1 - ^')y) {x* Jy*)] (31) 
^ - 1) .^^4 _ i)^*)(a;*ja;) + 2k\1 + K^)x*){y*Jx) + 2k\1 - K^)x){x*Jy*)] (32) 

KO 

where 6 is the scalar function 

6={1- K^f{x*Jx){y*Jy) + An^il + K^f{x* Jy){y* Jx) - - K^ f{x*Jy*){xJy) (33) 

At the moment we have fulfiled the reduction constraints (|22|) , ( ^3]) completely, but the con- 
straint eq.(0) has not yet been imposed. One can prove without loose of generality, that 
eq. (pl|) is verified, when the vectors x),y), fulfil 

y) = x). (34) 

Assembling eqs.(^) and eqs. (|3T|) -(^4|) all together one can find $ (eq.(p4|)) as a function of 

x). A, Aq 

$ = <l>(A,Ao,x)) (35) 

Now, we are interested in the value of this function at A = 0, since it gives a new solution 
g = $(0)g(o). One can check that 

g = g(o)-lil^[(l-«:4)(Fg(o) + g(o)F)(xVx)+ (36) 
2(1 + k^){k''F - g(o)Fg(o))(xVg(o)x) - 2(1 - K')Re{{K^H + g(o)i/g(o))(xVg(o)X*))], 

where 

F = Re(x*)(xJ), H = x){xJ 

Thus, we have parametrised the new solution g by a real number k and a complex vector 
x) of zero pseudolength. Since g satisfies the conditions (p!OD, it has the form of eq.(0) and the 
corresponding vector g), can be found by projecting eq.(^ on an arbitrary constant vector p). 
(For instance (1, 0, 0, 0)). We find in this way: 

,) = , . (37) 

^2{{pJp) - (pJgp)) 

Now, to construct a solution q) of the string equations (H, |^), one needs only to determine 
the dependence of x) and k on the variables ^,77. 
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3.3 Evolution of x) and in ^ and rj. 



Now the problem is to find the evolution of x) and k in ^ and rj. It follows from eqs.(|TT]), 
that 

+ Y3^<^'t/(o) = Y^A^*^' ^^^^ 

*. + Y^'^^{0) = Y^^*' (39) 

where f/, \^ are still undetermined functions of r], ^ which do not depend on A. Let us rewrite 
equations (^), (0) in the form 

t(-a5 + = (40) 

+ ^ (41) 

Consider the l.h.s. of (^Ol). It is a rational function of A with a pole at A = 1 and at A G 
{Ao, Ao^,Aq,Ao~^} , but the r.h.s. has only one pole at A = 1 . Thus to fit eq.(^OD we have 
to set the residues at Aq, Aq^, Aq and Aq"^ equal to zero. In fact, it is sufficient to require 
the vanishing of the residue at Aq only. All other residues will vanish due to the action of the 
reduction group. The condition 

resUo$(-9^ + ^)'^'"' = (42) 



yields 



and 



A(95--^)JA^ = (43) 

1 — An 



B A* B* U, 



+ + ^^^m - T^)JA' = (44) 



Ao — Aq Aq — Aq Ao — Aq 1 — a 







Both equations will be satisfied if 



{d^--^)x) = (45) 

i — Aq 



Thus, the simultaneous solution of eqs. (^3D and (0) is 

x) = ¥'\r],^-Xo)xo) (46) 
where Xq) is any complex constant vector of zero pseudolength, and Aq turns not to depend on 
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Moreover, the solution ( ^6]) of equation (|^) is also a solution of the equation 



Vi 



resUo$(-5, + Y^)*"' = 



(47) 



Finally, the 77,^ dependance of the vector x) is given by eq. f^Gf) . Together with eg .(^61) it 
gives the one soliton solution. 

The wave function '^{jj,^; A) corresponding to the one soliton solution (let us denote it by 
\l/i(?7,,^; A) can be regarded as a function of A, Aq, ^^^^^r],^, A), and xq) (see eqs.(^) and (p5|)). 

v|/i(ry, A) = $(A, Ao, J^(°Hr/, X^)Jx,))^^{r^, A) (48) 

The wave function corresponding to a n-soliton solution can be obtained recursively through 
the relation 

^^„(r/, e; A) = <I>(A, A„_i, J^n-i{r]. ^; A„_i) Ja;„_i))^„_i(r/, ^- A) , (49) 
and the corresponding solution of the chiral model g = \E'„(?7,,^;A = 0) will satisfy all the 



reduction conditions eqs.([l6|),([T7|) and (p!9D . 



4 The choice of the background solution 



Let us construct now explicit solutions by applying the above procedure . To begin with we 
shall consider a three dimensional de Sitter spacetime (D = 3). 

As a background starting solution we choose for simphcity the solution g(o) (a, r) found in 



ref. |T^. This solution corresponds to the trivial a = solution of the sinh-Gordon equation 
and it is given by. 



1 
71 



9(0)77 



1 
71 



/ sinhr \ 
coshr 

COSO" 

\ sino" / 

/ coshr \ 
sinhr 
-sino" 

\ COSO" j 



1 

V2 



ho) = 



1 

V2 



( coshr \ 
sinhr 
sincr 
\ — coscr y 

/ sinhr \ 
coshr 

— COSO" 

\ — sina ) 



(50) 



(51) 



For this solution, we have 



Qimi ~ 1{o)vn — "(0)5 Q'(o)5»7 — 0(0)5 — ^{o)v> "(o)»? — Qim- 
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(?(o)c^g(o),y) = -1, (g(o)-^g(o)) = 1, (6(o)J&(o)) = l, other (-J-) = O 



U{o) = 2g(o)^)(g(o) J - 2g(o)) (g(o)5 J, V(o) = 2g(o)^)(g(o) J - 2g(o))(g(o),,-'' 
where J is given by eq.(|). Let us define 

Q{^, V) = (?{o)^ g(o)r,, q{o),ko)) ' 
then we find by direct calculation that 



Q^JQ 



( 0, -1, 0, \ 
-1, 0, 0, 



0, 



0, 1, 



^ 0, 0, 0, 1 y 



U{o)Q = (0, 2g(o), 2g(o)g, 0), V(^o)Q = (2g(o), 0, 2q^o)v, 0) 



and 



/ 0, 0, 0, \ 
0, 0, -2, 
0, 2, 0, 

Vo, 0, 0, oy 

f 0, 0, 0, -1 \ 
0, 0, -1, 
0, 1, 0, 

V 1, 0, 0, J 

We have to solve two compatible equations for \1/ 



/O, 0, 
0, 0, 
2, 0, 

V 0, 0, 

/O, 0, - 

0, 0, 

1, 0, 

vo, 1, 



-1, 

0, 
0, 
0, 



-2, \ 
0, 
0, 

0, ; 

\ 

-1 





(0) 



^(0) .T,(0) 

1-A ' 



V< 



(0) 



1 + A 



Let us make the gauge transformation 



then 



i — A 

and r/) satisfies the following equations with constant coefficients 



i + A ' 
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( 0, 0, 

0, 0, 0, 

0, 0, 

V -1, 0, 0, 



\ 

-1 







/ 


0, 


0, 


0, 


-1 \ 






0, 


0, 







^ri — 






0, 


0, 







[ 


0, 


-1, 


0, 


) 



The solution will be defined on the Riemann surface T which covers twice the complex 



A: 



1 - A 



The points A 



1 + A 

±1 are the branching points of F. 



The fundamental set of solutions of eqs.(|56|)-(|58|) is given by: 



exp(/ir/ + /i 



-1 

—ifi 
1 

V 1 / 



-1 
1 

V 1 / 



It will be convenient to use the following linear combinations of the above solutions: 

S(/i,e,r7) = A(/i).n(/i,e,r7) 



where: 



A(/i) = diag(/i 1, 1) 







cosh 7, 


— sinh7, 


sin^, 


—cosO ^ 


1 




cosh 7, 


— sinh7, 


— sin6', 


cos6' 






sinh7, 


— cosh 7, 


— cos6', 


—sin 6' 




v 


— sinh7, 


cosh 7, 


— cos6', 


— sin6' / 
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and 

1 1 

9 = ^T] ^ , 7 = /ir^ + (65) 

/i /i 

Thus, ^, 77) as expressed by 

e, r/) = g(e, r7).A(^).n(/i, r^) (66) 
is the solution of the system (^) fulfiUing the constraint (|I^) for the solution g(o). That is : 

= 1,^,7^) = l-2\q^o) ><qio)\J. 



We want solutions periodic in a with period 2tc. We see from eq. (|62D -(|66D that we have 
hyperbolic functions on the argument 

^{fi-fi-') + ^{fi + fi-') (67) 

and trigonometric functions with the argument 

The solution to the a-periodicity condition requires to have 

/i = exp[za] (68) 

with real a and cos a and sin a to be rational numbers. The general solution is given by the 
Pythagorean numbers: 

— 2mn 

cos a = — 5- , sina = , m,n = integers (69) 

_l_ ji'i jji'i _|_ 

That is, 

m + in 

fx = , m, n = mtegers (70) 

m — in 

We get in this way solutions with period 27r(m^ + n"^) in a. Upon the rescaling: 

(T ^ (j(m^ + ra^) , T^T{m^ + n^), (71) 

we set the period to 2it. Notice that now the background solution g(o)(a, r) will be winded 
(m^ + n^) times around the origin in de Sitter space. 
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5 The soliton string solutions and their properties 



We have now all the elements to obtain the explicit expression for the solution \q{i],C,)) of 
the string eqs.(0)-(|). The explicit expression for the solution \l/'^°^(?7, ^; /iq) given by eq.(|66[) 
can be directly obtained by computing the indicated matrix multiphcation; Qij],^) is given 
by eqs.( ^0| - ^ ); A(/i) and n(/i,^,?7) are given by eqs. (|64D . This was done with the help of the 
computer program of symbolic calculation "Mathematica" . 

By projecting '^'^'^\r],^] /iq) thus obtained on a constant and complex null vextor |xo) , we 
have directly the vector The matrix g{C,,Ti) is obtained from eq.(PB|) also using "Mathe- 
matica". Finally, the explicit solution \q{ri,^)) is obtained by inserting g{C,, rj) in eq. (|37D . These 
string solutions depend on one complex parameter fi that depends on two integers n and m (see 
eq.([7y)), and one complex null vector xq) , that is, three complex independent numbers. Only 
two independent complex components remain in fact since g{C,,'ri) is homogeneous in xq). As 
can be seen in eq. (|36D , the change Xq) c Xq) where c is a complex number, leaves the solution 
invariant. The dependence on xq) is precisely like what happens for strings in D-dimensional 
Minkowski space-time, in which the solution depends on 2{D — 2) complex coefficients. They 
account for the {D — 2)-transverse degrees of freedom and for the two helicity modes (right 
and left movers) . Here, we are in three space-time dimensions and so we obtain two complex 
coefficients corresponding to the transverse degrees of freedom. 



It can be noticed, that the linear system (14) fulfilled by \l/(?7,^. A) is invariant under con- 



formal transformations on ^, 77 . Thus the dressing transformations do not generate conformal 
modes but only physical (transverse) modes. 

The vector xq) describes the polarization of the string; the integers m, n associated to the 
a periodicity, label the string modes. In Minkowski space-time, only one integer labels the 
right-modes and another one labels the left-modes. Here, we obtain two independent integers 
for each mode. Notice that our modes combine left and right movers in a non-linear and precise 
way. 

The resulting solution q{cr,T) = [q^ , , q"^ , q'^) is a complicated combination of trigonomet- 
ric functions of a and hyperbolic functions of r . From eqs.(|67|)-(f70|), we see that we have 
trigonometric functions on the arguments 

2mn — rt^ 

^ — 2 ^ — 2 

_|_ TfY^^ _|_ 

and hyperbolic functions on the arguments 

2mn - rn? — 

and r — - (72) 



That is, these string-solitonic solutions, do not oscillate in time. This is a typical feature 
of string unstability 0,0 which is present for strings in inflationary type backgrounds, i. e. 
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accelerated expanding like de Sitter, and in black holes. The new feature here is that the 
string does not oscillate in time, neither for r — nor for r ±00. It can be noticed that 
in decelerated expanding backgrounds, as it is the case in the standard Friedman-Robertson- 
Walker (FRW) expansion, string instability does not occur and the string behaviour is oscillating 
in r [0. This was recently confirmed for all values of r, in the FRW universe, where explicit 



string solutions has been found |T^. The non-oscillatory behavior in time can be understood 
from the fact that the string motion in de Sitter spacetime reduces to a sinh-Gordon equation 
with negative potential [Q. In D = 3, this is precisely the standard sinh-Gordon equation. 



whose potential unbounded from below (see fig. 1) is responsible of the instability. By defining 

e-('^'-) = -q^.q^ , (73) 

the string equations and string constraints (|^) in de Sitter spacetime can be reduced to the 
sinh-Gordon equation 

- aaa - + C"" = (74) 

Therefore, in order to find a solution in D = 3 de Sitter spacetime, one can start from a 
cr-periodic solution of the sinh-Gordon equation ( ffil) and insert it in the string equations (^: 

'dl -dl- e"('^'")] g(a, r) = (75) 

Then, one must solve the linear equation (^) in g((T, r) and impose the constraints (^ and 
(H). This is actually an alternative method to the dressing method, to obtain string solutions 
in de Sitter spacetime. 

The function e"^'^'^-' has a clear physical interpretation, as it determines the proper string 
size. The invariant interval between two points on the string, computed with the spacetime 
metric, is given by 

ds' = ^^dq.dq = ^ e-^"^'-) {da' - dr') (76) 
The energy density of the sinh-Gordon model here 

n = l[{^r + {^r]-2cosha{a,T) , (77) 

determines the potential 

V;// = -2cosha , (78) 

This potential has absolute minima at a = +oo and a = — oo. As the time r evolves, a{a,T) 
generically approach one of these infinite minima. The first minimun corresponds to an infinitely 
large string whereas the second one describes a collapsed configuration. That is, the string in 
de Sitter spacetime will tend generically either to inflate (when a —>■ +oo) or to collapse to a 
point (when a — > —oo). 

The background string solution g{o)(cr, r) given by eq.([50|) , corresponds to the sinh-Gordon 
solution a = 0. This means a string with finite constant proper size (equal to 1/H ). In the 
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sinh-Gordon model this corresponds to a particle at the maximun of the potential Vejf = —2 
and with zero velocity. 

Let us recall that for a given time qo, the de Sitter space is a sphere S"^ with radius 



R = jjJl + ql . For the background solution g(o)((T, r) given by eg. (^0]) , we have R{t) 



+ \ sinh^ r. As de Sitter universe expands for r — > oo , the string size e"*-'^'^-' = 1 remains 
here constant. This solution is probably unstable under small perturbations. 

It must be noticed that the integers (m, n) of the solitonic solutions (^) have the meaning 
of string winding. They label the different ways in which the string wind in the spatial compact 
dimensions (here 5*^). Notice that our string solutions do not oscillate in time in spite of the 
fact that we are in a lorentzian signature spacetime. (The dependence on r is hyperbolic). 

In figs. 3 and 4 we plot the one soliton solutions |g((T, r)) = {q^ , q^ , q^ , q^) found here. 
They show the three dimensional spatial projections {q^,q'^,q^) as a function of a, for a given 
polarization vector different values of m,n and two differents values of r. Figs. 3a and 3b 
show the same solution (n = 2,m = 1, x°) = (1, —1, .1, .li)) for two different values of r. Figs. 
4a and 4b show the evolution for a higher winding number {n = 5), and polarization vector 
xO) = (l,-l,l,^). 

For comparison, let us recall, that in D = 2 spacetime dimensions, in which the string 
motion reduces to the Liouville equation, the exact general solution is a string wound n times 
around the de Sitter space and evolving with it. The string covers n times de Sitter space which 
is here a circle S^. This solution is given by |]l4 



n 1 cos racT 9 smna , 

q" = - cot nr , g = , q= (79) 

sm riT sm nr 

< d < 27r , < r < 7r/n (80) 
The invariant interval between two points of the string 

ds^ = (da^-dr^) (81) 

sin nr 

exhibits the typical feature of string instability : in the asymptotic regions r —>■ 0^ and t ^ n/n, 
the proper string length blows up. We also see that the string does not have "enough time" to 
oscillate in one expansion time of the universe : the oscillation period of the string coincides 
with the expansion time of the universe. When the string accomplishes one oscillation, the 
universe has ended. 

In order to analyze the exact solutions |g(cr, r)) , it is convenient to use the coordinates 
{T, , X^) in this 2 + 1-dimensional de Sitter spacetime: 

gO = sinh HT +^ exp{HT) [ {X^f + {X^ ] (82) 

= cosh HT exp{HT) [ {X^Y + {X^ ] (83) 



18 



That is, 



= Hexp{HT)X^ , q^ = H exp{HT) , (84) 

-oo < T, X\X^ < +00. 



H H q^ + q^ H q^ + q^ 



The cosmic time T and the conformal time 77 are related by 

1 



T] = e-"^ ,-oo<r]<0 , (86) 



in terms of which, the line element takes the form 



ds' = -dT' + e'''' {dX'Y + idX'y =—— -{dr^Y + {dX^f + {dX'f . (87) 

L J H^Tj^ L -I 

We now analyze the properties and new features exhibited by these solutions. 

First of all, let us analyze the cosmic time coordinate T = T{(t, t). We have studied T as a 
function of r for different fixed values of a and viceversa. These functions have been obtained 
numerically for a wide family of solutions labeled by different values of the parameters n, m and 
\x'^) . We report here only two significative cases, which show the generic features, irrespective 
of the particular values of these parameters. 

Fig. 5 shows T as a function of r for the values of a indicated in the picture. We depict T 
for n = 4, m = 1 and a generic = (1 + i, .6 + Ai, .3 + .5i, .77 + .79i). 

In fig. 6 we depict T as a function of r for the solution with = 4, m = 1 and = 
(1, — 1). We see that our solution in the generic case describes actually five strings, as it 
can be seen from the fact that for a given value of T we find five different values of r .That is, 
r is a multivalued function of T for any fixed a . This is an entirely new feature for strings 
in curved spacetime. It has no analogy in flat spacetime where the time coordinate obeys 
[(9^ — 9^] T = , and therefore, using the conformal transformations 

cr±r = /±(cr'±r ) 

allows to choose the light-cone gauge, in which T (or a null-like combination of it) is proportional 
to r. In curved spacetime, this is not possible in general. Only in some geometries like 
shock- waves or gravitational plane waves, the light-cone choice is possible for all r [Q. In 
asymptotically flat spacetimes this choice is only possible asymptotically [||. When r is a 
univalued function of the time coordinate, the solution of the string equations describes only 
one string. This is the case in geometries where the the light-cone choice of gauge is possible for 
all r . In spacetimes, as de Sitter, where r is a multivalued function of the time coordinate, the 
solution of the string equations of motion and constraints describe a multi-string configuration. 
That is, each branch of r as a function of T corresponds to a different string. 
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This multiple number of strings arises as a consequence of the string dynamics in curved 
spacetimes, that is, from the coupling of the string with the spacetime geometry. Notice that 
here we have just free string equations of motion in curved spacetime. That is, interactions 
between the strings themselves, hke sphtting and merging, are not considered. We find that 
the geometry determines the simultaneous existence of several strings. They do not interact 
directly between them since they do not intersect. All the interaction is through the spacetime 
geometry. Notice that such phenomenom does not appear in D = 2, [eq.(|80D] where time is a 
monotonic and periodic function of r. This solution describes only one string in one period : 
< T < 71 /n. For others periods we get identical copies of the same string. This is not the case 
of the 2+1- dimensional solutions displayed in figs. 5-11. They describe five or three different 
strings. Five is the generic number of strings in our dressed solutions. This value five can be 
related to the fact that we are dressing a one-string solution (g(o)) with four poles. Each pole 
adds here an unstable string. 

Figs. 7 and 8 show the function r = r(cr, T) as a function of a -the different values of T are 
indicated in the pictures- for the above solutions (i.e. polarization vectors > and windings 
[n, m) the same as above). The function r = r((T, T) being periodic in a , it is plotted only for 
one period (27r) of a. We see that in addition to the period 2tt, another period in cr appears which 
depends on r. r = r((T, T) is a sinusoidal type function. It is more convoluted for small values 
of |r| in the neighbourhood of r = where several maxima and minima appear. As soon as the 
neighbourhood of r = is left, r = r((T, T) becomes very fast a regular sinusoidal-type function 
of a with a fixed period much smaller than 27r. (In all solutions studied here, r = T{a,T) 
reaches this asymptotic form for r ~ 5). The meaning of these small and large r behaviours 
will become more clear in connection with the evolution of the spatial coordinates and shape 
of the string. The small r behaviors are connected with the different (and complicated) ways 
in which the string winds at the begining of its evolution, while the r — oo uniform behavior 
is connected with the asymptotic configuration which is 'frozen' in comoving coordinates. The 
large r behaviour turns to be r-independent in comoving coordinates. 

Let us analyze now the spatial coordinates X^{a,T) and X^(cr, r) of this solution. Figs. 
9-10 show the time evolution of the three or five strings simultaneously described by this 
solution. In order to describe the real physical evolution, we eliminated r = r((T, T) from the 
solution and expressed X^(cr, r) = X^{a,T) and X^(cr, r) = X^(o", T) in terms of T. This was 
done numerically. Figs. 10 show the comoving coordinates (X^,X^) for different times HT. 
We see that for the fifth string, (X^, X^) collapse precisely as the inverse of the expansion factor 
e~^'^ , while the other four strings keep (X^,X^) constant in time (in Fig. 9, it is the third 
string that collapses). That is, the first string keep its proper size constant while the proper size 
of the other four strings expand like e^'^. These exact solutions display remarquably the string 
behaviour found asymptotically and approximately in refs.0. In summary, when (X^,X^) are 
smaller or equal than 1/H (the horizon radius), they contract to a point keeping the proper 
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amplitudes {e^^ X^,e^^ X"^) and proper size constant. When (X^,X^) are larger than 1/H, 
they become very fast constant in time, the proper size expanding with the universe itself as 
e^^ (string unstability). 

In terms of the sinh-Gordon description [see eqs.(|73|)-(|7^ and fig. 1] , this means that for 
strings outside the horizon, the sinh-Gordon function a{a, r) for most of the history is the same 
as the cosmic time T up to a function of a. We find combining eq.(5.14) in ref.[n] with eg. (pSf) : 



a{a, r) ^^=" 2HT{a, r) + log {2H^ [{A^aYY + (^^(a)')'] } + 0(e-^^^). (88) 

Here A^{a) and A'^{a) are the and X"^ coordinates outside the horizon. For T oo the 
string is at the absolute minimun a = +00 of the sinh-Gordon potential and possess an infinite 
size. 

The string inside the horizon corresponds to the maximun of the potential, a = 0. This 
is the stable string with contracting coordinates (X^,X^) and constant proper size, appearing 
in all the multi-string solutions found here. The value a = is the only in which the string 
can stay without being pushed down by the potential to ±00. This also explains why only one 
stable string appears : it is not possible to put more than one string at the maximun of the 
potential without falling down. The starting zero soliton solution a = we have dressed is a 
particular and very simple stable string. 

For degenerate choices of the number of strings reduces to three [see fig. 6]. For large 
positive T two of the strings (strings 1 and 2) are of the unstable type and one (string 3) is of 
stable type. In addition, strings 1 and 2 become identical in the infinite T limit. In fig. 9, we 
plot this solution for negative T. We see that string 2 is stable for T — * —00 (it has constant 
invariant size in such limit) whereas the invariant sizes of strings 1 and 3 collapse in this limit. 
In addition, there is an intermediate regime for \T\ < 3 where the comoving size of the strings 
decreases by a factor of about 10. 

The features above described are generically exhibited by our one-soliton multistring solu- 
tions independently of the particular initial state of the string. (Fixed by the values |x° > and 
(n, m)). 

It is interesting to see how the shape of the string becomes more symmetric for special values 
of >. For instance, a rosette shape or a circle with many festoons are clearly shown by figs. 
9-11. They correspond to >= (1,-1, i,l) with n = 4,6, respectively. These particularly 
symmetric vectors |x° > yield also degenerate solutions, in the sense that they contain only 
three different strings instead of five, as it happens in the generic case. 

We also see that for the symmetric inital conditions for the string state >, the function 
T = T{cr, T) becomes a perfectly symmetric periodic sinusoidal inside the period 27i, for all 
values of r (including small |r| ), and the additional very small period is practically the same 
for all |r| . 

The number of string windings and festoons is related to the frequencies in eq. (|69D and 
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express in terms of (n, m) . The a dependence is characterized by the frequencies [see eq. (|69D ] 

2mn m? — 

^ — 2~i — 2 ' = — 5— — 5- 

and the basic frequency 

1 2 
^0 = — ;5 o + odd , l^o = o + even . 

For m = 1, the highest available frequency is the sum ^1+^2 = " This highest frequency 
determines the small period in r = r(cr, T) as a function of a for fixed large r. That is, 

2.- "' + 1 



n2 + 2n - 1 
In addition, 

— — = + 2n — 1 for odd n, — — = [n + 2n — l)/2 for even n 

gives the number of festoons in the strings at a given T (see figs. 7-8 and 9-11). 

Strings propagating in de Sitter spacetime enjoy as conserved quantities those associated 
with the 0(3,1) rotations on the hyperboloid (^). They can be written as 

2iT I r2iT 







dcy{\q){q\~\q) {q\) J = - T da {U + V) 

I Jo 



In order to compute L it is convenient to relate U and V with f/(o) and V(o) using eqs.(|40|)-(^lD 
and the asymptotic behaviour of $(A) for A ^ 00 [see eq.( 



$(A) = 1 + ili^ + 0(-^) 



where C{r],C,) is a matrix. We then find 

U + V = U^o) + no) + 2C{7],0. 
Since C{ri,C,) is a periodic function of a, 

L = L(o) 



for all solutions considered here. We recall []15[ that only Lqi does not vanish for g(o), taking 
the value [see eq. ([7T|) ] : 

Lio = -Loi = (n^ + m'^)'K 
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Figure Captions: 

Figure 1: Effective potential corresponding to the sinh-Gordon model. 

Figure 2: Function HT{t, u), for fixed a — .41, for the n — 1 string solution in 1 + 1-de 
Sitter spacetime. 

Figure 3: Evolution of the string in (c, r) variables. The three projections {q^,q'^), (?^,?^) 
and (g^, q^) are shown for n = 2, m = 1 and r = 0.1 and 2 for >= (1, —1, .1, .li). 

Figure 4: Same as fig.3 for n = 5, r = 0.5, 3.85 and \x^ >— (1, —1, i, 1) 

Figure 5: Plot of the function HT{t), for two values of a, for n — A, |x° >= (1 + i, .6 + 
.4i, .3 + .5i, .77+ .79i). The function t{T) is multivalued, reveahng the presence of five strings. 

Figure 6: Same as fig. 5, for n = 4, >= (1, 1). Because of a degeneracy, there are 
now only three strings. 

Figure 7: r = t{(t,T) for fixed T for n = 4, \x^ >= (1, — l,i, 1). Three values of HT are 
displayed, corresponding to HT=0 (full line), 1 (dots), and 2 (dashed line). For each HT, three 
curves are plotted, which correspond to the three strings. They are ordered with r increasing. 

Figure 8: Same as fig. 7 for n = 4, x° >= (1 + i, .6 + Ai, .3 + .5i, .77 + .79i). a) The five 
curves corresponding to the five strings at HT=2. b) The five curves for three values of HT: 
HT=0 (full line), 1 (dots), and 2 (dashed line). 

Figure 9: Evolution of the three strings, for n = 4, >= (1, —1, i, 1). The comoving size 
of string (1) stays constant for HT < —3, then decreases around HT — 0, and stays constant 
again after HT — 1. The invariant size of string (2) is constant for negative HT, then grows 
as the expansion factor for HT > 1, and becomes identical to string (1). The string (3) has a 
constant comoving size for HT < —3, then collapses as e~^^ for positive HT. 

Figure 10: Evolution of three of the five strings for n = 4, |x° >= (1 + i, .6 + .4i, .3 + 
.5i,.77+.79i). 

Figure 11: Evolution of the three strings for the degenerate case n = 6, |a;° >= (1, —1, i, 1). 
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